Electronic Properties of Disordered Graphene Antidot Lattices 
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Regular nanoscale perforations in graphene (graphene antidot lattices, GAL) are known to lead 
to a gap in the energy spectrum, thereby paving a possible way towards many applications. This 
theoretical prediction relies on a perfect placement of identical perforations, a situation not likely 
to occur in the laboratory. Here, we present a systematic study of the effects of disorder in GALs. 
We consider both geometric and chemical disorder, and evaluate the density-of-states as well as 
the optical conductivity of disordered GALs. The theoretical method is based on an efficient algo- 
rithm for solving the time-dependent Schrodinger equation in a tight-binding representation of the 
graphene sheet [S. Yuan et al., Phys. Rev. B 82, 115448 (2010)], which allows us to consider GALs 
consisting of 6400 x 6400 carbon atoms. The central conclusion for all kinds of disorder is that 
the gaps found for pristine GALs do survive at a considerable amount of disorder, but disappear 
for very strong disorder. Geometric disorder is more detrimental to gap formation than chemical 
disorder. The optical conductivity shows a low-energy tail below the pristine GAL band gap due to 
disorder-introduced transitions. 

PACS numbers: 73.21. La,72.80.Vp,73.22.Pr 
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I. INTRODUCTION 



Pristine graphene has no band gap: the conduction 
and valence bands touch at the K and K'-points of the 
hexagonal Brillouin zone. This property, combined with 
the linear dispersion of the low-energy excitations leads 
to the spectacular electronic properties that graphene is 
so famous for— Nevertheless, the lack of a gap severely 
hampers many applications where a gap is needed to con- 
trol the flow of charges. This feature is further under- 
scored by the phenomenon of Klein tunneling: graphene 
carriers impinging on a potential barrier may experience 
refiectionless tunneling thus making their control even 
more difficult^. It is thus natural that many schemes 
have been proposed to create a gap in graphene: these- 
suggestions include etching extended graphene flakes into 
nanoribbons 5 -^, or by considering bilayer graphene in a 
transverse electric field 7 -^, or by using an external pe- 
riodic potential to modify the electronic properties so 
that a gap is formed. The external periodic potential 
may be caused by a number of agents, such as periodic 
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or strain^, or adsorption of adatoms in a regu- 



lar patter n i , or, as in this work, by a regular nanop- 
crforation of the pristine graphene sheet; this system will 
be referred to as graphene antidot lattices (GAL) 14 . 

The design principle^ behind the GAL was inspired 
by photonic crystals where pass and stop bands for 
light can be designed by drilling holes in the dielectric 
medium. GALs (and their constituents, single holes 
in graphene^ 5 -) have been studied theoretically with a 
large number of methods, ranging from a continuum 
description and tight-binding mcthod a 14 ' 16 to fully mi- 
croscopic DFT calculation s 17 ' 18 . Both electronic 



,19-22 



and thermal^ 3 - - — transport properties, as well as opti- 
cal propertie d 26 ' 27 have been discussed. Symmetry prin- 
ciples determining the existence or non-existence of the 
gap have been outline d 28 ' 29 . Most important, however, 
is the recent emergence of experimental techniques by 
which GALs can be fabricated. These fabrication meth- 
ods include, e.g., electron-beam etching 3 ^ - — , etch-masks 
based on self-assembled block co-polymers 3 ^—, nanoim- 
print technology 3 - 7 ., or nanoparticle depositio n 38 ' 39 . Most 
experimental papers have focused on the structural as- 
pects, but also a few transport experiments have been 
reporte d 30 ' 31 ' 33 ' 35 ' 39 . Indeed, transport gaps have been 
observed but so far they have been associated to dis- 
order induced localization instead of band-structure 
effect o 33 ' 35 . This highlights the importance of studying 
disorder in GALs: all fabricated structures contain disor- 
der, and one cannot (yet) control the exact geometry of 
the edges of the etched holes. It is thus vital to examine 
the robustness of the band-gaps against disorder, whether 
it be structural, geometrical or chemical. A study of this 
kind presents a serious computational challenge because 
the systems fabricated in the lab, where unit cells of the 
order of tens of nanometers can be achieved, arc com- 
putationally large involving tens of thousands of carbon 
atoms in the computational cell. Fully microscopic DFT- 
based methods cannot presently address such systems, 
and certain compromises must be made. 

In this paper, we perform a systematic study of the 
electronic properties of disordered GALs in the frame- 
work of a tight-binding model in a perforated honeycomb 
lattice of carbon atoms. We consider the most relevant 
kinds of disorder for these systems, namely a random de- 
viation of the periodicity and of the radii of the nanoholcs 
from the perfect array, as well as the effect of resonant 
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scatterers in the sample (like vacancies, adatoms, etc.) 
and the effect of non-correlated and correlated (Gaus- 
sian) on-site potentials. Within this scheme, the density 
of states (DOS) is obtained from a numerical solution of 
the time-dependent Schrodinger equation (TDSE) 4 ^, and 
the optical conductivity is calculated by using the Kubo 
formula for non-interacting electrons! 40 ' 41 

The paper is organized as follows. In Sec. |II]we present 
the details of the method. The effect of the different kinds 
of disorder on the DOS and the optical conductivity of a 
GAL is discussed in Sec. IIIII Finally, our main conclu- 
sions are summarized in Sec. IIVI 

II. MODEL AND METHOD 

We consider the following real space tight-binding 
Hamiltonian for a disordered GAL 

U = - ^2{Uja\bj + h.c) +y^ y Vic\ci, +H imp , (1) 

where a\ (hi) creates (annihilates) an electron on sublat- 
tice A (B) of the honeycomb graphene lattice, iy is the 
nearest neighbor hopping parameter and Vi is the on-site 
potential. In our model, the GAL is simulated by the cre- 
ation of an hexagonal array of circular holes of a given 
radius R, and a separation P = \f?>L between the cen- 
ters of two consecutive holes, where L is the side length 
of the hexagonal unit cellJ^ We thus label our GALs 
with the parameters {L, R}, in units of the graphene lat- 
tice constant a = \/3a ~ 2.46 A, where a ~ 1.42 A is 
the carbon-carbon distance. Another possible notation 
is [P, N], where N = \/3L — 2R is the neck width, defined 
as the smallest edge-to-edge distance between two neigh- 
bouring holes in the array^ Deviations of the GAL with 
respect to perfect periodicity are considered in our cal- 
culation in a twofold manner. First, we allow the center 
of the holes to float with respect to their position in the 
perfect periodic lattice (x,y) around (x ± lc, y ± lc) [see 
Fig. [T](a)]. Second, we let the radius of the holes to ran- 
domly shrink or widen within the range [R — rR, R + Tr] . 
as sketched in Fig. Q] (b). All along this paper, we will 
express lc and tr in units of a. 

The second term to the right of Eq. (flj accounts for 
a change in the on-site potential of the carbon atoms. 
A long-range potential for correlated impurities can be 
modelled with 

«. = f>»p(-^). w 

where N c is the number of the impurity centers, which 
are chosen randomly distributed on the carbon atoms, 
Vk is uniformly random in the range [— Vb,Vb] and d is 
interpreted as the effective potential radius. The value 
of iV c is characterized by the ratio n c = N c /N, where 
N is the total number of carbon atoms of the sample. 



A non-correlated short-range random potential can be 
obtained from the above equation with d — > 0, i.e., i>i is 
random and uniformly distributed, independently of each 
site i, in the range [— v r , +v r ]. The number of sites with 
nonzero pontential (N r ) is characterized as n r = N r /N . 

We further consider the effect of isolated vacancies in 
the sample, which can be regarded as an atom (lattice 
point) with an on-site energy t>, — ¥ oo or, alternatively, 
with its hopping amplitudes to other sites being zero. In 
the numerical simulation, the simplest way to implement 
a vacancy is to remove the atom at the vacancy site [see 
Fig. 01(c)]. 

If additional resonant impurities are present in the 
sample as, e. g., hydrogen adatoms, their effect is ac- 
counted for through the term Himp in Eq. (TTJ): 

i i 

where Sd is the on-site potential on the "hydrogen" im- 
purity (to be specific, we will use this terminology al- 
though more complicated chemical species can be con- 
sidered, such as various organic groups^) and V is the 
hopping between carbon and hydrogen atoms. 40 ' 42 ' 43 The 
spin degree of freedom, which contributes through a de- 
generacy factor 2, is omitted for simplicity in Eq. (JTJ) . All 
along this work, we fix the temperature to T = 300K. We 
use periodic boundary conditions in the calculations for 
both the optical conductivity and the density of states, 
and the size of the system is 6400 x 6400 atoms. 

Our numerical method is based on an efficient eval- 
uation of the time-evolution operator e -1 ^*, based on 
the Chebyshev polynomial representation^ - (In fact, any 
function of H can be evaluated with this method). We 
have thus access to the time-dependent state \<f(t)) = 
e -iUt |^ where \ip) is a random superposition of all the 
basis states in the real space, i.e. ) 40 ' 44 

M=5>cJ|0), (4) 



a-i are random complex numbers normalized as \ai\ = 
1, and |0) is the electron vacuum state. 

The numerical method has the advantage that an av- 
erage over different random initial states is not needed. 
This is because one initial state contains all the eigen- 
states in the whole spectru m 40 ' 44 . Furthermore, it is not 
necessary to average over different realizations of the dis- 
order, because the system contains millions of carbon 
atoms, and one specific disordered configuration contains 
a large number of different local configurations. As shown 
in Ref.[40], the results for different disorder configura- 
tions are essentially identical. 

Consider first the optical conductivity. We omit in our 
calculations the w = Drude contribution to the real 
part of the optical conductivity, so that the regular part 



Figure 1: Sketch of the different kinds of disorder considered, a) The center of the holes is shifted randomly with respect to 
the original position in the perfect periodic array (x,y) to a new position in the range (x ± lc,y ± lc) (lc = 2a). (Notice 
the different relative distance between the holes) b) The radius of the holes is randomly shrunk or enlarged within the range 
[R — tr, R+tr] (j-r = a). (Notice the different relative size of the holes) c) GAL with additional randomly distributed vacancies, 
signaled by the missing carbon atoms (n x = 1%). d) GAL with randomly distributed hydrogen adatoms, signaled by the red 
dots (rn = 1.75%). Notice that two other kinds of disorder are considered in the text, namely non-correlated and correlated 
long-range (Gaussian) changes in the on-site potentials, which are not sketched in this figure. 

sample area, / (%) = 1/ [e^^ - ^ + l] is the Fermi-Dirac 
distribution operator, and we use units such that h = 1. 
The time-dependent current operator in the a (= x or 
y) direction is J a (t) = e lHt J a e~ lHt . The Fermi-Dirac 
distribution operator / {%) is computed with the Cheby- 
shev polynomial representation, as mentioned above. As 
the next example, consider the overlap between the time- 



can be written a a 40 ' 45 
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e -0u _ 2 

Cq/3 (w) = lim / e~ Et am ujt 

x2Im {<p\f (H) J a (t) [1 - / (%)] Jp\<p) dt, 

(5) 

where j3 = 1/ksT is the inverse temperature, Q is the 
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Figure 2: DOS (top panels) and optical conductivity (bottom panels) for a {10, 6} GAL with geometrical disorder. In the first 
column we show (a) DOS and (c) <r(w) for a disordered GAL in which the center of the holes is shifted randomly with respect 
to the original position in the perfect periodic array within the range (a; ± lc, y ± lc), as sketched in Fig. Q](a). The different 
colors correspond to different values of lc (in units of a), as denoted in the inset of the figures. In the second column we show 
(b) DOS and (d) er(w) for a GAL where the radius of the holes is randomly shrunk or enlarged within the range [R — tr, R + tr], 
as sketched in Fig. Q](b). Different colors correspond to different values of ra (in units of a). 



evolved state \<p(t)) and the initial state \<p). The Fourier 
transform of this object yields the DOS of the system 

a = 40,44 

1 f°° 

P& = ^ J J st (vW{t))dt. (6) 

Finally, the quasi-eigenstate | $(£■)), which is a su- 
perposition of the degenerate cigenstates with the same 
cigenencrgy E, is obtained as the Fourier transform of 

\$(E)) = — dte* Et \ V {t)). (7) 

^ J-oo 

The quasi-eigenstate is not exactly an energy eigenstate, 
unless the corresponding eigenstate is not degenerate at 
energy E. However we can still use the real space distri- 
bution of the amplitude to examine the quasi-localization 
of the mode o 40 ' 41 ' 46 . Below we display several examples 
of all these objects. 



III. RESULTS AND DISCUSSION 

In this section we present the results and discuss the 
effect of the different kinds of disorder introduced in Sec. 
HT1 in the DOS and in the optical conductivity of GALs. 
As discussed in Sec. [TT] and sketched in Fig. [TJ we con- 
sider three main sources of disorder: geometrical disor- 
der, which is associated to deviations of the GAL from 
the perfect periodicity; resonant impurities, which can be 
associated to additional vacancies in the graphene lattice, 
or to adatoms deposited on the sample; and the effect of 
on-site potentials which can randomly vary within the 
sample. 

A. Geometrical disorder 

We start by considering the most generic source of dis- 
order in these kind of systems, which is the geometrical 
disorder. Uncontrollable fluctuations in the fabrication 
process lead to irregularities in the resulting antidot lat- 
tice, such as changes in the center-to-center distance of 
the etched holes, or in variations in the size of the holes. 
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Figure 3: Distribution of quasieigenstates of a {10, 6} GAL with geometrical disorder, where the radius of the holes is randomly 
shrunk or enlarged within the range [5.75,6.25]. The holes with unchanged radius (R — 6) are indicated by the red arrows. 
The quasi-eigenstates are calculated at energy E — 0.056t and E — 0.12t, corresponding to the states at the low energy peaks 
marked by the black and red arrows respectively in the DOS of Fig. b (tr = 0.25). The highest contribution to the quasi- 
eigensate at E = 0.056t (plot a) is concentrated around holes with radius R > 6, for which the first ring of atoms, as compared 
to perfect GAL with R = 6, has been removed. The quasi-eigenstate at E = Q.12t with highest amplitude (plot b) are localized 
around the holes with R — 6 (marked by the red arrows), as in perfect GAL. 



Examples of the geometrical disorder in the lattice are 
sketched in Fig. [T] (a) and (b), respectively. Let us con- 
sider first the effect of a random deviation of the relative 
distance among the holes on the DOS and cr(u)), as shown 
in Fig. (5Ja) and (c) respectively. In Fig. f^a) we see 
that, for the perfect periodic array (lc = 0, black line) 
a clear band gaps open up in the spectrum. Notice that 
the peaked structure of the DOS is due to the set of lo- 
cally flat bands which appear in the new band structure 
of the GAL as compared to the spectrum of standard 
grapheneJ^ If we now allow for a relative displacement 
among the nanoholes (lc 7^ 0), we observe that the gap 
shrinks but survives if lc is not too large, as seen by the 
red line of Fig. [2] (a), but eventually the gap closes for 
some critical value of lc, due to the lack of periodicity 
in the GAL, as it is the case shown by the green line of 
Fig. [5] (a) which correspond to lc = a. The disorder 
affects the optical conductivity in the following manner. 
As it can be seen in Fig. [5] (c), for the perfect periodic 
case (lc = 0) cr(uS) = up to ui = A, where A w 0.2t 
(for the case considered here) is the gap opened due to 
the antidot array. Because A decreases as we increase 
lc, the threshold for optical transitions is reduced and 
for lc = 0.5a we observe a finite optical conductivity for 
oj > 0.05t. Finally a(uj) > at any frequency for an even 
larger amount of disorder as, e. g., lc = a (green line), 
for which the gap of the GAL has completely collapsed. 

A similar effect on the electronic properties is observed 
if instead of randomly changing the relative separation 
between the antidots, their size is varied within some 
range, as sketched in Fig. Q] (b). The results of our sim- 



ulations for this kind of disorder are shown in Fig. [2] (b) 
and (d) for the DOS and optical conductivity, respec- 
tively. One observes that the DOS presents an increas- 
ing number of peaks as r# is increased. These peaks are 
associated to states with a large amplitude circling the 
antidots, and their energy depends on the radius of the 
antidot. For the pristine GAL (tr = 0, black line) all 
antidots have the same radius which leads to the peaks 
at E/t ~ ±0.12 [signaled by a black arrow in Fig. [5] 
(b)]. The finite width of the peak is due to a coupling 
between the antidots, and the weak splitting reflects the 
van Hove singularities at the edges of these quasi-one 
dimensional bands. Peaks at higher energies originate 
from states that are not tightly localized around the an- 
tidots, but have a larger amplitude all over the sample. 
If the radius of the antidots is varied we observe that, 
apart from the peak discussed above shown by the black 
arrow, part of the spectral weight is transferred to new 
peaks that correspond to localized states at different en- 
ergies, around antidots of different radii. Some examples 
are shown by the red and green arrows in Fig. [5] (b). 
This behavior is illustrated by the spatial distribution of 
the quasi-eigenstates shown in Fig. [3] There we show, 
for the case of = 0.25, a small section of the lattice 
studied in our simulations, with a real space distribution 
of the amplitude of the quasi-eigenstates |$(i?)| 2 . Fig. 
[3] (a) shows that, for the energy E/t rj 0.056 [marked 
by the red arrow in Fig. [5] (b)], the large amplitude of 
the states is around antidots with R > 6, for which the 
innermost ring carbon atoms of the antidot has been re- 
moved. However, in Fig. [3] (b) we see that at an energy 
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corresponding to the first mode of the undistortcd lattice 
[E/t ps 0.12, shown by the black arrow in Fig. [2] (b)] the 
states are localized at the edges of holes with a radius 
corresponding to perfect GAL. Those antidots are shown 
by red arrows in Fig. [3j 

The effects of this kind of disorder on the optical con- 
ductivity are shown in Fig. [2] (d). As r# is increased, 
we observe optical processes of lower and lower energy 
contributing to o~(ui), due to optical transitions between 
the localized states around holes of different sizes. This 
suggests that photoluminescence spectroscopy can be an 
useful tool for the characterization of the GALs. 



B. Resonant impurities 

The next source of disorder that we consider is the 
effect of resonant scattcrers. Resonant impurities can 
be understood as vacancy atoms in the sample, or as 
hydrogen or other organic molecules (CH3, C2H5, etc) 
adsorbates which bind to a single carbon atom, chang- 
ing its hybridization from sp 2 to sp 3 — A sketch of 
a GAL with a certain amount of vacancies or hydrogen 
adatoms randomly distributed is shown in Fig. [T](c) and 
(d) respectively. The main effect of the resonant impuri- 
ties in graphenc membranes is the creation of " mid-gap" 
states at the Dirac point . 41 ' 42 Therefore, if some amount 
of these kind of impurities is present in the GAL, a zero 
energy flat impurity band is expected to appear in the 
middle of the gap. This is indeed what we obtain in 
our calculations, as it can be seen by the E ps peak 
in the DOS plots shown in Fig. 2J As we discussed in 
Sec. |nl this kind of disorder is accounted for in our cal- 
culations through the term "Hi mp in Eq. (|3|), with the 
band parameters V ps 2i and e c i ps — 1/16, as obtained 
from ab initio density-functional theory.— The DOS of 
a GAL with different amounts of vacancies and hydro- 
gen adatoms, randomly distributed, is shown in Fig. [4] 
(a) and (b) respectively. We observe that, apart from a 
slight deviation from the Dirac point (E = 0) of the posi- 
tion of the hydrogen adatoms impurity band (due to the 
finite value of the energy Ed), as compared to the E = 
energy of the mid-gap band due to vacancies, the effect 
of these two kind of defects in the spectrum is very simi- 
lar. In the two cases, the quasi-localization of the newly 
created states leads to an almost flat band which does 
not affect the rest of the energy spectrum away from the 
Dirac point (apart from some smearing of the peaks in 
the DOS). 

As a consequence, the main contribution to the optical 
conductivity is obtained, as in the clean limit, for inter- 
band processes with an energy ui ps A, as it can be seen 
in Fig. S] (c)-(d) . However, due to the transfer of spectral 
weight to the mid-gap states, there is some finite a(u) for 
energies smaller than the threshold defined by the energy 
gap A, with an appreciable peak at w ps A/2 ps O.lt. 
This contribution is due to the new optical transitions 
from the impurity band to the conduction band, which 



are activated for u> > A/2. 

C. Short- and long-range potential disorder 

Another kind of disorder that can be considered is a 
shift of the on-site potentials at a given lattice points, 
which can lead to a local shift of the chemical poten- 
tial. This contribution is accounted for by means of the 
second term of the Hamiltonian (TT]). This kind of disor- 
der can be of extraordinary importance. For example, if 
the atoms in sublattices A and B have opposite strength 
of the on-site potential v r , then a gap of size A = 2v r 
is opened in the spectrum^ Here we consider, depend- 
ing on how the defects are distributed over the lattice 
sites, a correlated or a non-correlated disorder. In the 
case of a short-range and non-correlated potential disor- 
der, the nonzero on-site potentials are taken to be uni- 
formly randomly distributed over the sample within a 
range [—v r ,v r ]. The results for the DOS and the optical 
conductivity with v r = Zt are shown in Fig. [5] (a) and 
(c) , respectively. We observe a broadening of the peaks 
in the DOS, accompanied by a transfer of spectral weight 
to the gaped regions. 

Next, consider the long-range correlated disorder given 
by Eq. In standard graphene, this kind of disorder 
leads to regions of the graphene membrane where the 
Dirac point is locally shifted to the electron (V4 < 0) or 
to the hole (14- > 0) side with the same probability. This 
leads to some finite DOS at zero energy. Our calculations 
for GALs in the presence of a Gaussian potential disorder 
are shown in Fig. [5](b) and (d), for the DOS and optical 
conductivity respectively. One observes similar qualita- 
tive effects in the spectra as compared to the short-range 
non-correlated random potentials. In particular, there is 
a small but appreciable contribution to the optical con- 
ductivity at low frequencies, due to the transfer of states 
to the gapped region. 



IV. CONCLUSIONS 

In conclusion, wc have presented a systematic study of 
the effect of disorder in GALs. We have used a tight- 
binding model in a perforated honeycomb lattice of car- 
bon atoms. The DOS has been calculated from a numer- 
ical solution of the TDSE, whereas the optical conduc- 
tivity has been obtained by using the Kubo formula for 
non-interacting electrons. We have considered the most 
generic sources of disorder in these kind of samples: ge- 
ometrical disorder such as random deviation of the peri- 
odicity and of the radii of the nanoholes from the perfect 
array, as well as the effect of resonant scatterers in the 
sample (e.g., vacancies, adatoms, etc.) and the effect of 
non-correlated and correlated (Gaussian) on-site poten- 
tials. In order to have a qualitative understanding of the 
effect of the different kinds of disorder on the samples, 
wc have applied the method to one representative case, 
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Figure 4: DOS (top panels) and optical conductivity (bottom panels) for a {10,6} GAL with resonant impurities. In the first 
column we show (a) DOS and (c) a{ui) for a GAL with a random distribution of vacancies, as sketched in Fig. [1] (c). The 
different colors correspond to different amounts of missing dangling bonds, as denoted in the inset of the figures. In the second 
column we show (b) DOS and (d) u{oj) for a GAL with hydrogen adatoms, as sketched in Fig. [l](d). Different colors correspond 
to different percentage of adatoms in the sample. 



namely a {10,6} GAL. However, we emphasize that this 
scheme is completely general and applicable to any set of 
parameters {L,R}. 

Our results show that the gap is rather robust against 
geometrical disorder, and only a large deviation of the 
antidot array from the perfect periodicity leads to a nar- 
rowing and eventually closing of the energy gap. We 
obtain localized states encircling the antidots, the en- 
ergy of which depends on the radius of the hole. The 
presence of additional resonant scattcrers, as vacancies 
or adatoms, leads to the creation of midgap states. The 
existence of this impurity band is reflected in the opti- 
cal conductivity, which now extends to energies smaller 
than the gap energy A, due to disorder activated opti- 
cal transitions from the impurity band to the conduction 
band. However, the main contribution to a(u>) still cor- 
responds to transitions with an energy of the order of 
A. Finally, the presence of non-correlated or of Gaussian 
potential disorder leads to a smearing of the peaks in the 
DOS, as well as to the transfer of spectral weight to the 
gapped region. Contrary to the effect of resonant scat- 
terers, the presence of potential disorder does not create 



a zero-energy band with a prominent peak in the DOS 
at E = 0, but instead a DOS that grows smoothly as a 
function of energy within the gapped region. As a conse- 
quence, the optical conductivity also grows slowly from 
until it reaches its maximum contribution at the energy 
of the gap. Therefore, photoluminescense spectroscopy 
experiments could be useful for the characterization of 
the GALs. 
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Figure 5: DOS (top panels) and optical conductivity (bottom panels) for a {10,6} GAL with on-site potential disorder. In 
the first column we show (a) DOS and (c) <t(u)) for a GAL with a non-correlated random distribution of short-range potential, 
which can take the values within the range [—v r ,v r ]. The different colors correspond to different concentrations of disorder, as 
denoted in the inset of the figures. In the second column we show (b) DOS and (d) er(aj) for a GAL with a long-range Gaussian 
potential disorder. The potential is given by Eq. ((2|, where Vk is uniformly random in the range between — Vo and Vo, and d 
is the effective potential radius. (See text). 
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